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This  paper  presents  the  system  description  and  organization  of  maps. 
: Map  Assisted  Photo  interpretation  System.  MAI’S  is  a large 
integrated  database  system  containing  high  resolution  aerial 
photographs  digitized  maps  and  other  cartographic  prirducts. 
combined  with  dclailcrKII)  descriptions  of  man-made  and  natural 
features  in  the  Washington  l).C.  area.  A classification  of  image 
database  systems  into  three  models  is  also  presented.  I hese  models  are 
the  Image  Database  (ID)  model,  lire  Map  Picture  Database  (MI’D) 
model  and  die  Image/Map  Database  (imu)  model.* 


1 . Introduction 

This  paper  presents  the  system  description  and  organization  of  MAPS, 
the  Map  Assisted  Photo  interpretation  System,  maps  is  a large 
integrated  database  system  containing  high  resolution  aerial 
photographs,  digitized  maps  and  other  cartographic  products, 
combined  with  detailed  3D  descriptions  of  man-made  and  natural 


2.  Background 

Our  early  motivation  for  investigating  image  databases  was  as  a 
component  of  a complete  image  understanding  system.  We  had  only  a 
vague  idea  of  what  capabilities  it  should  have,  hut  we  thought  that  it 
should  represent  "idealized  segmentations”  of  an  image,  where  die 
labeling  of  die  segments  was  in  fact  the  "scene  interpretation".  It 
should  relate,  or  compare  machine  generated  segmentations  to  dris 
model,  and  provide  die  user  with  a qualitative  and  quantitative 
performance  measure  of  die  machine  segmentation.  We  attempted  diis 
with  die  MIDAS  system1  ■ using  the  segmentation  results  for  a set  of 
Pittsburgh  city  scenes  generated  by  the  argos'-  4 system  Ihc  results  of 
the  performance  analysis  of  the  scene  segmentation  were  less  than 
encouraging.  While  we  could  give  quantitative  analysis  of  the 
segmentation  and  labeling  by  the  ARGOS  system,  the  qualitative  results 
were  couched  in  the  original  (subjective)  hand  segmentations.  It  was 
difficult  to  qualitatively  distinguish  between  alternative  machine 
segmentations,  since  the  relative  importance  tor  cost  function)  of 
missing  or  mislabeled  regions  or  broken  boundaries  for  different 


features  in  the  Washington  D.  C.  area. 

This  paper  discusses  three  major  topics.  First,  a classification  of 

dilTcrcnl  models  of  database  systems  for  cartographic  applications  is 

presented  together  with  a discussion  of  their  inherent  strengths  and 

limitations.  I hese  models  arc  the  Image  I kinbase  (id)  model,  lire  Map 

Picture  Database  (mhd)  luodel  and  the  Imagc/Map  Database  (imd) 
* ' • * 

model.  Second  -wc  argue  lor  Ihc  tuility  of  (he  lin.igc/M.ip  Database 
model,  discuss^lasks  and  present^!  general  description  of  the  model. 


regions  was  not  represented  in  the  segmentation.  How  to  perform  such 
an  evaluation  is  still  an  open  research  problem.  Also,  although  we  had 
a database  of  18  high  resolution  color  images  of  Pittsburgh,  we  had  no 
general  mechanism  to  relate  one  to  another,  except  through  analysis  of 
the  hand  segmentations  and  the  names  given  to  buildings,  roads,  rivers, 
and  other  features  in  the  scene.  However,  in  the  process  of 
implementing  and  using  MIDAS  we  did  learn  a great  deal  about  image 
database  organization  and  symbolic  representation  of  scene 
descriptions. 


I bis  model  describes  components,  facilities  and  techniques  that  should 
be  present  in  such  a sy  stem,  and  a range  of  tasks  that  can  he  supported 
by  the  model.  Finally,  we  descrilnxthc  maps  system  in  terms  of  our 
(IMD)  model,  and  discuss  three  applications  which  utilize  and  integrate 
image,  terrain,  and  map  data  in  a powerful  manner.  Wc  also  discuss 
what  wc  leave  learned  during  the  implementation  of  the  MAI’S  system, 
some  ideas  on  the  proper  interfaces  between  components,  where 
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modularity  should  he  ;ichicvcd.  and  point  to  future  work. 


Dm  research  was  sponsored  by  ihc  Defense  Advanced  Research  Projects  Agency 
(DOD).  ARPA  Order  \o  3597  monitored  by  the  Air  1‘orcc  Asiorncs  laboratory  Under 
Contract  I 33M5-R1-K-I539  Ihc  vww>  and  conclusions  contained  in  this  document  are 
those  of  the  author  and  should  not  he  interpreted  as  representing  the  official  policies, 
either  expressed  or  implied,  of  the  Defense  Advanced  Research  Projects  Agency  or  the 
US  Government 


Wc  decided  to  look  .it  map-guided  image  interpretation  and  began  to 
assemble  an  aerial  photograph  database  of  the  Washington.  I).  C.  area. 
Using  this  imagery,  we  felt,  we  could  quickly  generate  a map  database 
that  would  allow  us  to  explore  image  analysis  of  complex  aerial 
photographs  using  a simple  map  database  that  constrained  where  to 
look,  and  what  to  look  for.  This  idea  of  map-guided  segmentation  was 
not  new.  Ihc  llAWkl  Yi  system5  and  succeeding  "road  cxpcrt"b  7 were 
based  on  similar  ideas,  and  use  of  world  knowledge  had  been  a well 
accepted  paradigm  in  image  interpretation.  However,  we  wanted  to 
focus  tin  more  general  capabilities,  to  represent  large  scale  spatial 
organizations  normally  encountered  in  complex  urban  sc  ncs.  ITic 
generation  of  the  map  database  turned  out  to  he  a much  harder 
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problem  than  wc  initially  estimated,  and  it  quickly  became  the  focus  or 
our  research.  In  retrospect.  I believe,  it  was  exactly  the  right  problem  to 
work  on.  and  although  there  is  still  much  to  do  in  the  area  or 
iinnge/map  dalnhnses,  wc  oow  have  the  right  tools  and  understanding 
to  begin  to  tackle  the  original  problem.  This  work  has  direct 
application  in  three  areas: 

• pholo-intcrprctalion:  representation  of  world  knowledge 
for  image  understanding. 

• situation  assessment:  a spatial  expert  for  decision  support 
systems. 

• cartography:  toward  digital  map  generation  and  use. 

3.  Classification  of  Databases 

There  has  been,  over  the  last  ten  years,  a perceived  need  for 
organizing  and  structuring  image  and  map  data  for  cartographic 
applications.  It  has  been  difficult  to  compare  various  capabilities  and 
limitations  of  systems  because  there  weie  few  common  denominator, 
by  which  systems  could  be  compared.  Systems  reported  in  the 
literature  could  loosely  be  categorized  either  as  research  vehicles,  or 
production-oriented  systems  for  particular  well  defined  subtasks  of  the 
general  cartographic  problem8' 9' 10.  Research  vehicles  generally  had  a 
high  degree  of  organizational  complexity  tested  on  very  small  scale 
databases.  Systems  used  in  production  environments  tended  toward 
simple  niodcis  running  very  large  stale  dauibasts.  f uilliCr,  Ktii!c  fix 
tasks  being  performed  involved  the  analysis  of  aerial  or  satellite  data,  it 
is  often  unclear  whether  the  image  tlaut  was  an  integral  part  of  the 
resulting  database,  or  simply  used  for  data  acquisition.  One  example  is 
the  development  of  digital  filing  systems  that  store  facts  about  a laigc 
iiuuihci  of  Im.igcs  without  the  actual  image  J:«t  Uve  hr-d 

example  of  such  a system  is  the  l-.ROS  Data  Center  database 
ri’nUWllMfl  by  dfe  l"S  lVp<  of  ih<*  Imorior  This  database  has 
approximately  2xl06  frames  of  l.andsat  imagery  and  5x10^  frames  of 
aircraft  (aerial  mapping)  photography.  Users  may  specify  an  area  of 
interest  by  geodetic  point  or  rectangular  area  and  stib-sclcct  those 
frames  based  on  time  of  year,  cloud  cover,  type  of  sensor  and  a a 
scene  quality  rating.  However,  the  actual  frames  of  data  are  stored  on 
high  density  magnetic  tape.  Similar  situations  exist  in  map  producing 
organizations  such  as  the  United  Stales  Geological  Survey  (L'SGS)  and 
die  Defense  Mapping  Acency  DMA. 

One  notable  exception  is  described  in  Kondo  curl.11  where  an  image 
database  using  l.andsat  imagery  was  integrated  with  map  descriptions 
for  geographic,  natural,  and  cultural  features.  Features  can  be 
displayed  superimposed  on  the  image  data,  and  imagery  could  be 
indexed  by  geodetic  location  or  by  feature  name.  There  arc  limitations 
such  as:  the  image-to-map  correspondence  was  based  on  a fixed 

decomposition  ofiandsal  data  into  a laliludc/longitudc  grid  at  a map 
scale  of  1:50000:  the  spatial  relationships  between  features  were  entered 
manually;  and  the  overall  complexity  of  the  image  and  map  database 
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was  small.  Nevertheless,  this  represents  an  ambitions  new  direction  for 
die  development  of  land-use  systems  using  l.andsat  imagery. 

In  this  discussion  of  database  systems  for  cartographic  and  situation 
assessment  applications,  wc  are  assuming  that  die  following  minimal 
capabilities  hold:  (I)  on-line  display  of  digital  imagery  and  map  data, 
and  (2)  ability  to  query  interactively  about  attributes  of  the  imagery  and 
map.  The  following  is  onr  classification  of  the  capabilities  or  three 
models  which  wc  can  use  to  compare  various  existing  systems  or 
approaches.  These  models  are  the  Image  Database  (111)  model,  die  Map 
Picture  Database  (MI’O)  model  and  die  Image/Map  Database  (IMD) 
model. 

3.1 , Imago  Databases 

flic  Image  Database  model  (ID)  is  the  simplicsi  and  most  common 
database  model.  It  is  organized  to  relate  attributes  about  the  sensed 
image  such  as  sensor-type,  acquisition,  cloud  cover,  or  geodetic 
coverage’**.  These  databases  generally  do  not  represent  die  content  or 
thu  scene,  hut  rather  attributes  or  the  scene.  When  the  semantics  of  the 
scene  arc  present,  the  location  of  cartographic  features  are  represented 
in  the  image  (pixel)  coordinate  system.  'I  bis  poses  obvious  limitations 
to  tlie  application  of  relevant  knowledge  from  other  images  or  from 
external  sources,  since  there  is  no  general  mechanism  to  relate  map 
feature  position  between  images  that  overlap  in  coverage  or  to  an 
external  map.  Although  the  features  represented  may  appear  to  be 
man-oriented,  is  is  difficult  to  compute  general  gcomcuic  properties 
using  the  image  raster  as  the  coordinate  system. 

Although  relational  database  techniques  have  been  applied  to  die  in 
model,  wc  feel  these  techniques  arc  not  appropriate  to  spatial  database 
organizations  for  several  reasons.  First,  using  the  basic  <aitribute ; 
nilutO  tuple  to  represent  vector  Iims  of  unip  coordinate  data  requires 
that  all  of  the  primary  key  attributes  be  duplicated  in  cacli  relation, 
aioec  dime  o lib  fntjll.iflimi  tlWHitilg-  nduod  lkts 

order  pairs)  as  a primitive  attribute  in  a , elation.  Further,  the  relational 
database  operations  such  as  union,  intersection,  join,  project,  arc  not 
good  primitives  fiir  implementation  or  inherently  geometric  operations 
such  as  containment,  adjacency,  intersection  and  closest  point. 
Operations  such  as  fcatuie  intersection  arc  reduced  to  searching  for  line 
segments  which  share  the  same  pixel  position.  Finally,  in  any  large 
system,  a logical  partitioning  of  the  database  must  be  performed  in 
order  to  avoid  extensive  and  often  unnecessary  search  when  performing 
spatial  operations.  Partitioning  is  difficult  to  achieve  itl  relational 
systems  since  the  relational  model  restricts  itself  to  homogeneous  (only 
one  record  type)  sequential  sets.  Previous  work  advocating  such 
organizations  did  not  address  die  issues  of  system  scale,  and  focused 
more  on  issues  of  query  languages  using  relational  models  for 
geographic  databases  than  the  actual  construction  of  complex 
systems12'1314.  When  measured  by  die  number  of  images,  image- 
based  features,  and  by  the  complexity  of  die  relationships  represented, 
these  systems  were  quite  simplistic. 

’’’using  night  annotation  such  as  the  center  point  and  corner  points  not  using  general 
imagc*to*niap  correspondence 
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3.2.  Map  Picture  Databases 

I lie  Map  IVuire  I )aiabasc  model  (mid)  describes  databases  that  arc 
generated  h\  digili/ing  cartographic  products,  such  as  pre  existing  maps 
and  charts  I licse  databases  arc  attractive  in  environ  rents  where  paper 
maps  have  played  a large  role  in  planning  .ind  analysis.  Ilierc  arc. 
however,  some  major  limitations  to  spatial  systems  based  on  digitized 
eat  tographic  products.  first.  in  the  original  map  production,  spatial 
ambiguity  lias  been  rectified  by  the  cartographer  in  a manner  that  is  not 
often  rcvcrsablc.  I be  cartographic  process  involves  simplication 
(generalization ).  classification  (abstraction),  and  symbolization  of  real- 
world  ambiguity.  Constraints  imposed  by  the  scale  of  die  map  often 
determine  which  woild  features  can  be  depicted  despite  the  desirability 
of  portraying  a complete  spatial  representation.  Therefore,  map  icon 
and  sy  mbology  placement  may  not  be  as  accurate  as  the  original  source 
material.  Since  the  deduction  of  the  actual  spatial  arrangement  of 
objects  from  an  iconic  representation  is  an  open  problem.  Ml’IJ's 
represent  chaos  masquerading  as  rationalized  order.  Hie  key  issue  is 
that  vim's  are  pictures  of  a map  (however  detailed)  rather  Ilian  the 
underlying  map  structure  and  spatial  organization.  Although  the 
graphics  display  of  Mi’ll  appears  to  convey  a great  deal  of  semantic 
information,  that  impression  is  a result  of  die  human  observer,  not  a 
reflection  of  an  underly  ing  map  representation. 

When  a map  is  digitized  into  a map  picture,  another  subtle 
simplification  recurs.  The  digitization  process  results  in  a map  image 
on  a rectangular  gi  id  whose  size  is  generally  limited  eidier  by  custom  or 
as  ..n  artifact  of  die  digitization  process.  Common  limitations  arc 
scanner  resolution,  maximum  size  of  image  raster,  and  the  physical  si/c 
of  source  map.  One  popular  representation  is  to  subdivide  regions  of 
the  nap  picture  into  a regular  decomposition  such  as  quad-tree15  l6.  or 
k-d  tree17.  Ihe  implementation  of  this  representation  is  greatly 
simplified  in  Mi’ll  models  since  one  no  longer  has  to  contend  with 
positional  ambiguity  of  map  features  because  of  the  cartographic 
process  outlined  above,  and  the  discrete  nature  of  die  digitization 
process. 

One  common  use  for  die  MID  model  is  in  geographic  information 
systems  for  kind  use  and  urban  planning.  In  diese  systems,  aggregrate 
values  such  as  population  of  an  area  and  crop  yield  of  an  area  are 
computed.  I he  scale  of  the  original  map  becomes  the  limiting  f.ictor  for 
accuracy  in  information  computation.  However,  the  grain  of 
computation  is  usually  large  enough  that  these  inaccuracies  arc  not  a 
practical  problem.  Incremental  update  of  the  database  due  to  new 
residential  ami  industrial  areas  and  the  coiH.-omil.int  loss  of  rural  areas  is 
a difficult  problem  since  database  update  requires  careful  map  editing 
tools  not  usually  associated  with  these  Mi’ll  systems. 

A recent  trend  has  been  to  take  existing  Mi’ll  databases  and  add  a 
map  feature  database  component  usually  relational  to  describe 
attributes  of  various  features.  We  believe  that  augmenting  traditional 
Mi’ll  databases  with  semantic  information  has  merit  in  those 
enviroments  where  analysis  is  being  performed  by  humans,  since 


information  synthesis  is  not  a requirement  of  the  database  system. 
However,  once  such  a system  is  in  place,  there  is  a tendency  to  attempt 
to  automate  analysis  functions  requiring  spatial  interpretation,  and  the 
generation  method  of  the  MI’li  model  has  several  drawbacks  for  use  in 
photo-interpretation,  situation  assessment,  and  cartography.  I be  chief 
problems  arc  the  method  of  generation  as  outlined  above,  the  lack  of 
semantic  information  about  map  features,  and  the  requirement  that  a 
map  exist  at  the  appropriate  level  of  detail  for  the  area  under 
consideration.  Ihe  IMI1  model  discussed  in  the  following  section 
addresses  these  issues. 

3.3.  Image/Map  Databases 

Hie  Image/Map  Database  model  (IMIl)  relates  map  features  to 
image  database  through  camera  models.  It  therefore  has  the  capability 
to  describe  relationships  between  features  acquired  from  different 
images  through  the  map  database.  This  capability  is  in  contrast  to  the 
image  database  model  where  the  feature  descriptions  can  only  be 
related  if  the  descriptions  come  from  the  same  image. 

Since  the  map  database  is  built  directly  from  aerial  imagery  in  the 
IMIl  model,  the  resolution  / accuracy  issue  is  a function  of  the  ground 
resolution  of  the  imagery  , the  intrinsic  position  measurement  error  due 
to  camera  model,  ground  control,  etc.  rather  than  an  artifact  of  the  map 
depiction  scale  as  in  the  MPD  model.  A greater  variety  of  feature 
descriptions  is  possible  since  they  are  not  restricted  to  those  that  can  be 
portrayed  in  a cartographic  product,  f urther,  the  complexity  of  a 
particul.it  feature  description  is  independent  of  any  particular  task 
requirement  and  can  represent  a rich  set  of  attributes,  semantic 
interpretations,  and  knowledge  from  diverse  sources.  I his  flexibility  is 
a key  clement  for  map  data  representation  as  we  look  toward  spatial 
database  systems  with  applications  in  cartographic  production,  expert 
photo- interpretation,  and  situation  assessment. 

However,  just  as  the  cartographer  must  resolve  ambiguity,  so  the 
spatial  database  must  be  able  to  represent  inconsistency  in  a consistent 
manner,  l-or  example,  errors  in  correspondence  between  images  and 
the  geodetic  model  cause  the  same  point  on  the  earth  to  be  given  a 
different  geodetic  position,  ie.  when  viewed  from  different  images  the 
same  geodetic  point  produces  a different  world  position.  If  this  point  is 
on  a common  boundary  between  two  features,  say  a political  boundary, 
there  should  be  ambiguity  as  to  which  region  the  point  is  in.  By  tile 
same  token,  if  two  large  residental  areas  are  found  to  intersect  because 
of  positional  uncertainly,  and  the  result  of  the  intersection  is  several 
small  polygonal  areas,  the  IMIl  model  should  be  able  to  rectify  this 
ambiguity.  Iliis  rectification  might  take  the  form  of  a symbolic 
relationship  that  indicates  that  the  residential  area  share  a common 
boundary,  while  maintaining  the  ability  to  represent  the  original 
errorful  signal  data.  Since  the  original  data  is  maintained  in  the 
database,  the  symbolic  relationships  do  not  have  to  be  static,  for 
example,  these  relationships  can  be  dependant  on  attributes  sinuliar  to 
those  used  by  cartographers  when  they  perform  simplification  and 
generalization.  Ihe  link  from  the  symbolic  interpretation  back  to  the 
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original  source  data  is  not  possible  in  Ml’lJ  systems. 

3.3.1.  Spatial  Knowledge 

I lie  imii  model  goes  us  the  tools  to  construct  our  map  database  from 
"first  principles"  and  tie  together  partial  spatial  knowledge  at  different 
lescls  of  detail.  Ibis  is  possible  because  individual  map  features  may  be 
specified  directly  from  source  imagery.  I bis  capability  is  precluded  by 
the  derivative  nature  of  die  Ml’l)  model.  That  is.  it  is  difficult  to 
assimilate  new  and  possibly  errorful  knowledge  because  of  the 
mismatch  between  die  new  errorful  data  and  the  cartographic 
rectification  of  ainbigious  data. 

I he  representation  of  a multiple  levels  of  detail  paradigm  is  often 
invoked  as  a pari  of  a coarse-fine  or  hierarchical  matching  stiaiegy  in 
Ullage  processing  and  interpretation.  Coven  the  scale  and  digitized 
ground  resolution  of  an  image,  the  IMI)  model  can  generate  a map 
description  that  will  suppress  any  features  that  would  be  urn  small  to  be 
recognized,  with  remaining  descriptions  at  the  appropriate  level  of 
detail.  I bis  technique  is  more  than  camera  scaling  and  transformation, 
smee  the  criterion  for  "loo  small"  can  be  an  attribute  of  the  map  feature 
itself  Consider  die  map  feature  description  of  a university  campus.  At 
some  level  of  detail  corresponding  to  pixel  ground  resolution  distance 
ii. Kin.  features  such  as  playing  fields,  dormitories,  instructional 
buildings  and  offices,  access  roads,  and  campus  greenery  arc  ^ ad 
individually  distinguished.  I 'sing  spectral  properties  of  the  features 
and  spatial  relationships  between  these  features,  we  can  determine 
those  feature  boundaries  that  arc  likely  to  be  muddled,  and  those  with 
sufficient  detail  to  be  recognized. 

Ibe  multiple  level  of  detail  paradigm  need  not  be  applied  in  a 
homogeneous  manner.  For  example,  tasks  such  as  decision  aids  for 
photo  intelligence  may  require  high  resolution  detail  to  support 
analysis,  but  low  resolution  detail  to  establish  overall  context.  A large 
scale  spatial  organization  containing  urban,  residential,  and  rural  areas 
will  require  flexibility  to  represent  the  high  feature  density  and 
complexity  in  the  urban  area  as  well  as  significantly  lower  density  in 
rural  areas. 

Flexible  knowledge  acquisition  is  necessary  because  in  photo- 
interpretation.  situation  assessment,  and  cartography,  world  knowledge 
is  inherently  fragmented.  Knowledge  fragmentation  in  these  domains 
.irises  from: 

• methods  of  know  ledge  acquisition 

llicre  arc  diverse  sources  of  kntmlcgc  th*ti  arc  used  to 
acquire  map  feature  information.  Some  of  the  most 
common  are  direct  measurement  from  imagery,  old  maps 
and  charts,  sketches,  and  collateral  data. 

• task  requirements 

If  the  task  requirement  is  to  support  radar  scene  simulation. 


then  elevated  roads  are  significant,  and  road  networks  in 
general  are  not  significant.  If  the  task  is  to  support  map 
generation  at  a particular  scale  (say  1:50000).  the  feature 
size  density  may  determine  whether  it  is  directly  portrayed, 
generalized,  or  omitted  entirely.  There  arc.  of  course,  well 
defined  rules  that  govern  these  decisions,  but  they  arc 
generally  not  consistent  across  a w idc  range  of  map  scales. 

• specialization  in  feature  extraction 

There  is  a certain  amount  of  specialization  in  cartographic 
and  situation  assessment  activities.  Analysts  may  specialize 
in  a particular  area  of  the  world,  be  know  ledgablc  in 
hydrology,  geology,  local  construction  customs,  or  political 
matters.  In  the  production  of  large  scale  maps  it  is  rare  to 
find  map  generalists,  although  tins  may  not  be  true  for  low 
level  feature  extraction  activities.  Ibis  specialization  tends 
to  fragment  knowledge,  and  is  often  given  as  a justification 
for  building  database  systems  that  provide  access  to  a wide 
range  of  map  knowledge  and  may  have  general  capabilities 
for  know  ledge  sy  nthesis. 

Ibe  IMI)  model  methodology  provides  a mechanism  for  feature 
unification  in  a cohesive  framework.  It  provides  a framework  to  relate 
symbolic  descriptions  to  their  original  data  sources.  It  is  not  tied  to  a 
particular  cartographic  representation  nor  to  limitations  of  cartographic 
production. 

4.  The  Database  Problem  in  Image  Interpretation 

I he  database  problem  has  been  addressed  in  a variety  of  ways  in 
systems  that  perform  image  analysis  and  interpretation.  However,  it 
has  rarely  been  pursued  as  a separate  research  problem.  One 
explanation  for  this  is  that  portions  of  general  database  represcnlion  arc 
often  embedded  in  tire  experimental  image  processing  systems  and 
become  highly  tuned  to  the  application.  Ibis  is  sometimes  a result  of 
system  performance  issues,  or  case  of  task-specific  implementations, 
but  often  it  is  a result  of  not  recognizing  llie  database  problem  as  a 
separate  issue. 

It  is  difficult  to  give  a precise  analy  sis  of  the  use  of  map  databases  in 
image  interpretation,  since  the  detailed  organizations  of  experimental 
systems  are  rarely  available.  However,  there  arc  several  recent 
examples.  Work  at  SRI  used  a map  database  of  road  intersections  to 

construct  a camera  model  in  the  iiawkfyi  and  subsequent  "road 

..  S h 7 

expert  systems 

Ibe  ARCjOS*' 4 system  used  a digitized  city  plan  map  and  elevations 
for  buildings  to  build  a .11)  graphics  model  of  downtown  l*ittsburgh. 
Ibis  model  was  directly  compiled  into  a knowledge  network 
representation  which  described  size,  shape  and  relative  positions  of 
buildings,  roads,  rivers,  and  bridges  for  an  arbitrary  view  point. 
Although  it  was  not  lied  to  a geodetic  grid,  it  was  a general  map  model. 


’“'for  example*  roads  preserve  linear  propen.es  until  Die  G«l>  approximately  equals 
ihc  width  of  the  road 
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Recent  work  .It  Hughes"'  based  on  the  ACRONYM  system  developed 
by  llrooks  and  Hitifordw  uses  image  registration  to  a geographic  model. 
Ibc  system  uses  pre-selected  regions  of  interest  and  attempts  to  locate 
and  identify  pre-defined  object  instances  within  dicse  areas. 

ACRONY  M is  currently  the  best  example  of  a model-based  system 
that  incorporates  viewpoint-insensitive  mechanisms  in  terms  of  its 
model  description.  Its  recognition  process  is  to  map  edge-based  image 
properties  to  instances  of  object  models.  In  the  domain  of  aerial  photo 
interpretation,  results  have  been  reported  for  the  recognition  of  a small 
number  of  models  (.1)  for  wide-bodied  jets  in  aerial  photographs.  It  is 
not  clear  how  map  knowledge  would  be  directly  integrated  into  die 
acronym  framework,  but  one  could  speculate  that  it  could  be  added 
by  a method  similar  to  the  work  at  Hughes  described  above. 

Matsuyama20  21  has  demonstrated  a system  for  segmentation  and 
interrelation  of  color-infrared  aerial  photographs  containing  roads, 
rivers,  forests,  and  residential  and  agricultural  areas.  It  uses  rules  to 
make  assignments  based  on  region  adjacency  and  multi-spectral 
properties.  Iltesc  rules  nuke  use  of  informal  map  knowledge  but  do 
not  directly  use  a particular  map  to  guide  interpretation.  It  generates 
good  descriptions  of  a variety  of  fairly  complex  aerial  scenes  getting  a 
great  deal  of  constraint  from  the  multi-spcctr.nl  data. 

In  his  recent  thesis.  Selfridge22  proposed  using  adaptive  threshold 
selection  for  region  extraction  by  hislogramming  and  region  growing 
using  an  image-based  "appearance  model”.  Although  die  work 
describes  feature  positions  and  shapes  in  terms  of  pixel  descriptions,  it 
is  not  difficult  to  imagine  a more  general  map-based  approach  dial 
would  result  in  the  automatic  generation  of  constaints  to  his  adaptive 
operators. 

At  C\1l . Herman  ! has  demonstrated  the  feasibility  of  incremental 
acquisition  of  31)  scene  descriptions  from  stereo-pair  aerial 
photographs  in  the  MAI’S  database  in  die  31)  Mosiac  project.  Ibis 
system  requires  a known  stereo  camera  model  but  uses  no  a-priori 
knowledge  about  the  scene  other  than  weak  geometric  assumptions 
about  urban  env  ironments. 

5.  The  Image/Map  Database  Model 

In  dus  section  we  discuss  four  classes  of  la  ’b  it  are  common  to 
photo  interpretation,  situation  assessment  and  cartography.  We  then 
list  some  criteria  by  which  one  can  evaluate  the  strengths  and 
limitations  of  database  systems.  These  criteria  arc  not  exhaustive, 
rather  they  point  to  four  areas  that  should  be  present  in  IMD 
implementations  and  system  capabilities  in  each  of  the  areas. 

5.1.  Tasks  lor  Image/Map  Database 

In  this  section  we  give  a classification  of  tasks  that  arc  common  to 
applications  in  photo-interpretation,  situation  assessment  and  digital 
cartography  systems.  The  four  tasks  ate  selection  of  image,  terrain,  or 


map  data  based  on  attributes  of  die  data,  spatial  compulation  of  map 
feature  relationships,  semantic  computation  of  map  features,  and 
synthesis  of  imagery,  terrain  and  map  data. 

1.  Selection 

The  selection  task  requires  that  the  IMD  system  he  able  to 
select  from  a potentially  large  set  of  database  entities  based 
on  attributes  of  image,  terrain,  and  map  database  features. 

The  selection  task  does  not  require  imagc-to-inap 
correspondence,  and  is  the  task  normally  performed  by  it) 
model  systems.  Tor  example: 

• select  imagery  with  particular  intrinsic  characteristics: 
sensor,  scale,  date,  cloud  cover,  processing  history 

• select  map  features  based  on  symbolic  description, 
partially  specified  description,  similarities  in  image 
acquisition 

2.  Spatial  Computation 

Spatial  computation  is  ubiquitous  in  cartographic,  situation 
assessment  and  photo-interpretation  tasks.  An  l\tl>  system 
must  provide  tools  to  compute  common  spatial 
relationships  such  as  containment,  closest  point,  adjacency, 
and  intersection.  One  issue  is  how  to  structure  the 
environment  in  order  to  constrain  search  and  thereby  avoid 
unnecessary  compulation.  Consider  four  views  of  the  same 
problem: 

• given  a geodetic  area,  which  images  cover,  or  partially 
cover  this  area 

• which  roads  can  be  found  within  the  image 

• which  images  contain  this  building 

• given  an  image,  find  all  images  which  overlap  it 

3.  Semantic  1 amputation 

There  are  a number  of  tusks  that  require  more  than  basic 
spatial  computation,  or  where  the  appropriate  spatial 
operation  depends  on  the  meaning  of  the  map  objects.  Arc 
there  intrinsic  high-level  properties  of  map  features  that  we 
can  extract  from  basic  spatial  geometry  that  give  a meaning 
to  the  feature?  Semantic  computation  needs  to  be 
investigated  as  we  develop  more  complex  spatial  databases. 

Tor  example,  what  is  the  semantics  of  intersection'  for  the 
follow  mg  pairs  of  map  objects? 

• intersection  of  two  roads 

• intersection  of  bridge  and  river  description 

• intersection  of  a building  and  a road 

4.  Synthesis 

One  goal  of  any  database  system  should  be  to  bring 
together  diverse  sources  of  knowledge  into  a common 
framework.  Sy  nthesis  is  the  generation  of  new  information 
using  a new  method  of  presentation,  computation,  or 
analysts.  Tor  example: 
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• cartographic  superposition  of  inap  data  on  newly 
acquired  image 

• 31)  display  of  terrain  and  cultural  features  from  map 
database  including  man-made  structures,  political 
boundaries,  neighborhoods,  arbitrary  collections  of 
physically  realized  features 

• to  predict  spatial  (location)  and  structural 
(appearance)  constraints:  where  to  look  and  what  to 
look  for  based  of  task  knowledge,  previous 
experience,  or  expectations 

• a spatial  framework  within  which  to  embedd  task- 
specific  knowledge 

5.2.  Criteria  for  Image/Map  Database 

In  this  section  we  list  some  criteria  that  can  be  uses  to  evaluate 
database  systems  in  four  general  areas.  These  areas  arc  imagc-lo-map 
correspondence,  map  feature  representation  spatial  computation,  and 
database  synthesis. 

1.  Image-to-Map  Correspondence 

• can  the  it  relate  image-based  features  to  a map 
coordinate  system 

• can  these  features  be  projected  onto  new  imagery 
using  the  correspondence  mechanism 

• chat  capabilities  exist  for  incrementally  updating 
feature  descriptions  based  on  updates  to  the  camera 
model,  or  to  intrinsic  changes  to  the  feature  self. 

3.  Representation 

• what  are  the  capabilities  for  feature  representation: 
what  complex  spatial  relationships  can  represent:  how 
is  inconsistency  recognized  and  handled 

• can  the  user  describe  features  and  associated 
attributes  in  a flexible  manner:  what  is  the  variety  of 
attributes. 

• can  the  representation  accommodate  map-based 
information  coming  from  a variety  non-imagery 
sources 

• what  is  the  relationship  between  the  representation  of 
signal  and  symbolic  data 

• what  synthesis  tasks  docs  the  representation  support 

3.  Spatial  Computation 

• docs  the  system  support  dynamic  spatial  queries 

• what  spatial  relationships  docs  the  system  compute 
directly  from  the  underlying  data,  which  relationships 
arc  specified  by  the  user,  how  do  they  interact,  how 
docs  one  maintain  consistency 

• what  mechanisms  arc  available  to  partition  the  search 
space  when  computing  spatial  relationships 


4.  Database  Synthesis 

• imagery,  terrain  and  map  data  arc  components,  each 
with  an  appropriate  representation,  operation 
semantics,  and  utility:  in  what  ways  does  the  database 
support  synthesis  of  these  components 

• what  concrete  lacks  requiring  synthesis  arc  performed 

6.  MAPS  Overview 

In  the  previous  sections  we  have  attempted  to  raise  issues  of 
Image/Map  Database  organization,  tasks  and  capabilities.  In  this 
section  we  will  discuss  the  MAPS  system  components  capabilities.  We 
will  only  briefly  describe  those  aspects  that  have  been  reported  on  in 
other  papers.  Our  latest  work  in  the  area  of  hierarchical  organization, 
decomposition,  and  search  is  reported  beginning  in  Section  6.6.  New 
work  in  map  feature  semantics  is  discussed  in  Section  6.7.  Tor  a more 
detailed  description  of  the  image  segmentation  program  (Section  6.1.2) 
and  the  image-to-map  correspondence  program  (Section  6.3)  see 
McKcown24.  For  a detailed  description  of  the  CONCI-TTMAP  database 
see  McKcown25.  Appendix  I contains  a nearly  complete  list  of  the 
programs  associated  w'ith  each  system  component. 

6.1 . BROWSE:  Interactive  Image/Mop  Display 

liROws I 2tl  is  an  interactive  window -based  image  display  system.  It 
provides  a common  interface  to  all  of  the  MAI'S  system  components  to 
display  results  of  queries,  graphical  prompts  for  interactive  image-to- 
map  correspondence,  superimpostion  of  map  data  on  imagery,  and 
other  similar  functions.  While  often  viewed  as  an  application  issue,  a 
flexible,  functional  user  interface  is  critical  for  building  more  complex 
tools,  iikovvsi  provides  the  user  with  a window -oriented  interface, 
which  greatly  increases  the  effective  spatial  resolution  of  the  frame- 
buffer. and  provides  multiple  processing  contexts  which  allow  users  to 
manipulate  dynamically  the  size,  level  of  detail,  and  visibility  of 
imagery. 

6.1.1.  Window-based  Display 

We  have  applied  and  extended  the  bit-map  window  paradigm  to 
handle  high  resolution,  multi-bit  per  pixel  digitized  images.  However, 
due  to  neat  < . an  order  of  magnitude  difference  in  the  amount  ol  data 
needed  to  perform  screen  updates  and  due  to  processing  limitations 
found  in  most  frame-buffer  architectures,  many  of  the  solutions  used 
for  single  bit  per  pixel  displays211  are  not  suitable  for  direct 
implementation.  A detailed  discussion  of  the  design  and  organization 
of  the  w indow  manager  appears  in  McKcown  & IX-nlingcr26. 

Besides  the  display  of  imagery,  we  have  found  the  window 
representation  to  be  useful  as  a communication  mechanism  between 
MAI’S  components,  to  invoke  image  processing  programs,  and  to 
retrieve  and  display  the  results  of  such  processing.  All  maps 
components  (sec  Appendix  I)  that  display  imagery,  map  data  or 
graphics  use  the  littowsi  window  mechanism  for  display  and 
communication.  For  example,  die  interactive  image  correspondence 
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program.  cokki  s.  uses  die  window  mechanism  lo  automatically  display 
landmark  image  fragments  and  lo  create  a high  resolution  window 
containing  the  approximate  position  of  the  landmark  ground  control 
point  to  cue  die  user,  fit tac  contains  a collection  of  image  processing 
routines  that  can  he  imoked  on  liKOWSl.  windows  simply  by  specifying 
the  window  name.  liKOWSl  routines  use  the  window  name  to 
determine  the  image  name,  resolution,  and  rectangular  image  bounds. 

I his  information,  along  with  parameters  specific  to  the  particular 
processing  operation,  are  passed  to  the  image  processing  routine.  The 
results  of  the  operation  can  he  displayed  in  a new  window. 

6.1.2.  Interactive  Image  Segmentation 

sruMl  vi  is  an  interactive  image  segmentation  program  which  uses 
the  liROWSI  window  facility  to  provide  an  interface  to  our  frame  buffer. 
Users  can  extract  image-based  descriptions  of  map  features,  edit 
existing  features,  and  assign  symbolic  names  to  the  features.  S1GMI  S I 
produces  a standard  formal  [Si  GJ  file  that  is  used  throughout  the  MAPS 
database  to  represent  image-based  descriptions  of  point,  line,  and 
polygon  geometric  data.  Database  routines  discussed  in  Section  6.S  arc 
available  to  convert  the  JS1G]  description  lo  a map-based  description 
IDJ]. 

6.2.  Image  Database 

Hie  MAPS  system  currently  contains  approximately  100  digili/cd 
images,  most  of  which  are  low  altitude  aerial  mapping  photographs. 
Iypic.il  ground  resolution  distances  (iiRli)  arc  120cm*.  160cm‘.  and 
600cm  per  pixel,  lire  imagery  is  mainly  comprised  of  three  data  sets 
taken  in  1974.  1976  and  1982.  In  addition  to  aerial  mapping 
photographs,  we  have  several  digitized  maps  including  a USGS 
topographic  map.  and  tour  guide  maps,  l-igurc  1 gives  the  current 
status  of  the  MAI’S  Washington  D C.  image  database.  Although  we 
have  several  I ands.il.  Skylah  and  high  altitude  aerial  photographs  taken 
over  the  Washington  D.C.  area,  we  have  focused  our  work  on  those 
images  that  prov  ide  the  greatest  ground  detail. 


CLASS 

NUMBER 

IMAGf 

SCALt 

DATABASE 

RASTER 

COMMENTS 

ASC • 74 

25 

1 36000 

2048.2048.8 

Aerial  napping 

BW 

WGl  76 

37 

1 12000 

2200*2200  » 8 

Atrial  napping 

BW 

AER  79 

? 

1 124000 

2288.2288.8 

Color  infrared 

ASC  8? 

29 

1 60000 

2300*2300  *8 

Aerial  n.ipping 

BW 

MAP  71 

1 

1 24000 

4096*4096*8 

USGS  topo  nap 

MAP  74 

1 

1 160030* 

4096.3880i8 

D.C.  region  nap 

MAP  79 

1 

l 16000* 

4096*4096*8 

Tourist  guide  i 

nap 

* not  car  toqrapli veal  ly  accurate. 


figure  I:  MAPS:  Image  Ihilabasc  Component 

6.2.1.  Generic  Image  to  File  Mapping 

Ibe  MAPS  system  uses  a gcncnc  naming  convention  to  refer  to 
images  in  the  database.  Ibe  generic  name  is  a unique  identifier 
assigned  to  the  image  when  it  is  integrated  into  the  database.  For 
example,  HC58M7.  DCHTO  arc  representative  generic  names  that 


correspond  to  flight  line  annotation  on  the  photographic  film.  All  types 
of  image  access  that  require  the  filesystem  name  of  die  image,  or 
require  associated  image  database  files,  use  the  generic  name 
mechanism  to  construct  the  appropriate  physical  file  name.  It  is 
possible  lo  change  the  logical  and/or  physical  location  of  imagery  by 
updating  the  generic  name  file  or  to  add  another  image  to  die  database. 
As  we  move  to  larger  image/map  systems  this  naming  isolation  allows 
us  to  construct  a database  that  can  be  distributed  over  multiple 
Ibe  decoupling  of  name  with  physical  or  logical  location  fits  well  with 
name  server  organizations  usually  employed  with  such  distributed 
systems. 

I lie  following  table  lists  the  database  files  associated  with  each  active 
image  in  the  maps  database.  Kadi  is  accessible  using  the  generic  image 
name. 

• |t , I si  ki<  ] imagc-to-filc  system  mapping 

- contains  the  file  system  location  of  the  database  image 

- identifies  which  reduced  resolution  images  arc  computed 
and  available  for  hierarchical  display 

• [SHI  1 scene  description  file 

- contains  image  specific  information:  source,  dale,  time  of 
day.  raster  size,  digitization,  image  scale,  geodetic  comer 
points,  camera  information 

• |t  Ol  ] iinage-to-map  coefficients  file 

- contains  camera  model  coefficients,  error  model, 
polynomial  orders  solved,  best  correspondence  (default 
polynomial  order) 

• independent  coefficients  for  <latiludc>.  <longitudc>. 
<iinage  row  > <image  column> 

• |t  OK)  correspondence  pairs  file 

- mapping  of  ground  control  points  to  image  point 
specification 

- lists  of  landmark  names  and  their  geodetic  position 
combined  with  image  pixel  position  of  landmark  specified 
by  user 

• [I IV I'J  hypothesized  landmark  file 

- lists  of  landmark  names  which  are  within  the  image 
geodetic  coverage,  but  were  not  used  to  perform  image-map 
correspondence 

6.2.2.  Image-Based  Segmentations 

MAPS  maintains  several  types  of  image  segmentations  and  map 
overlay  descnptions  associated  with  each  image  in  the  database.  These 
scgmentaions  either  are  feature  descriptions  generated  using  the  image 
as  the  base  coordinate  system,  or  the  projection  of  map  features  onto 
the  image  using  map-to-image  correspondence,  or  segmentations  from 
other  images  registered  to  the  image.  In  the  latter  ease,  image  to-map 
correspondence  is  used  to  register  the  two  images.  Users  can  point  to 
segmentation  overlay  features  using  the  display  interface  in  BROWSF 
and  com  I pim  xp.  identify  the  segmentation  feature  name  and  retrieve 
ns  image  and  geodetic  coordinates.  For  the  [tit  Mssuj)  and 
|<  omi  pish  segmentation  descriptions,  the  name  of  the  segmentation 
feature  is  used  to  retrieve  the  .ism Mated  in  All  (see  Section  6.4)  or 
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concttimap  description.  The  following  table  is  a list  of  imnge 

segmentations  associated  with  each  image  in  the  datahase. 
Segmentations  that  require  msip  correspondence  for  their  generation 
can  he  atitoniatically  rccieatcd  when  image  camera  model  is  updated. 

• IIIANHSPG]  hand  (human)  segmentation 

- collection  of  all  hand  segmentations  performed  on  this 
image 

• |llCO.MPSPG|  composite  hand  segmental  ion 

- collection  of  all  features  in  the  |Hasiisi-G]  database  that 
arc  spatially  contained  in  this  image 

• iMACMSl-.tj]  machine  segmentation 

• collection  of  all  machine  segmentations  performed  using 
the  image 

• [MCOMI’SP.G]  composite  machine  segmentation 

- collection  of  all  features  in  the  |MACIISIG]  database  that 
are  spatially  contained  in  the  image 

• |iii..vissi:g]  dims  map  overlay 

- all  features  from  the  di  ms  digital  feature  analysis  database 
dial  are  spatially  contained  in  the  image 

• ICONCMI’ISI  G]  CONCT-ptmap  map  overlay 

- all  features  front  the  conciptmap  database  that  arc 
spatially  contained  in  the  image 

• icovi  usi  Gl  image  coverage  overlay 

■ all  images  whose  area  of  coverage  is  overlapped  or  wholly 
contained  within  die  image 

6.3.  Image-to-Map  Correspondence 

The  MAI’S  system  uses  an  interactive  image-to-map  correspondence 
procedure  to  place  new  imagery  into  correspondence  with  the  map 
database.  It  has  three  major  components:  a landmark  database,  a 
landmark  creation  and  editing  program,  and  an  interactive 
correspondence  program.  The  process  of  landmark  selection, 
description,  and  interactive  correspondence  has  been  described  in  detail 
in  McKcown24. 

6.3.1.  Landmark  Database 

MAPS  maintains  a database  of  approximately  200  geodetic  ground 
control  points  in  the  Washington  D.C.  area.  Landmarks  arc  acquired 
using  USGS  topographic  maps,  bill  in  principle  can  lie  integrated  from 
any  source  that  provides  accurate  geodetic  position 

<liililmlc/biigiliilc/clemliwi>.  Users  can  query  the  database  to  lind 
landmarks  hv  name,  within  a geodetic  area,  or  the  closest  landmark  to  a 
geodetic  point.  Landmark  features  arc  also  integrated  into  die 
CONtTTIMAP  database  and  can  lie  found  using  die  <ralc-derimlion> 
attribute  (sec  Section  6.5.2)  of  a concept  role  schema. 

6.3.2.  LANDMARK 

I andmakk  is  an  interactive  tool  used  to  generate  new  landmarks, 
iCfl1  'll  Hi  ■fiTjnwrill*  The  following 

information  is  maintained  by  l andmakk  to  support  landmark  database 
access. 


• |l  dm]  landmark  name  directory 

■ associates  the  list  of  landmark  names  with  their  gcodclic 
position 

- sorted  for  spatial  proximity 

- partial  name  matching  also  provided 

• (m  Y|  landmark  lext  description 

■ contains  a detailed  text  description  of  tile  location  of  die 
landmark  and  general  factual  properties  of  die  landmark 

- stores  die  location  and  name  of  the  associated  image 
fragment  fdc  |l  IMG),  and  replicates  the  gcodclic  position 
from  klin  file 

• |l  IMG]  landmark  image  fragment 

• contains  a high-resolution  image  fiagmcni  which  clearly 
shows  the  ground  control  plant  and  scene  context  around 
die  point 

6.3,3.  CORRES 

corres  is  an  interactive  image-io-map  correspondence  program.  It 
uses  the  I1ROWSI-:  window  interface,  die  landmark  database,  and 
image  database  routines  to  interactively  liiiiid  an  image-tn-map 
correspondence.  Once  an  initial  guess  of  the  corner  points  is  performed 
and  the  |coit)  and  (COl!|  files  have  been  created  in  die  image  database, 
CORRi-s  aniomatically  suggests  new  possible  landmark  points  using  the 
image  database  [HYP]  Hies.  The  landmark  database  [i.lMG]  files  arc 
used  lo  display  the  ground  control  point  when  the  user  selects  it  from 
the  lisi  of  hypothesized  points. 

6.4.  dlms:  An  External  Database 

The  ability  ui  rendezvous  with  externally  generated  map  databases  is 
a key  capability  in  order  to  integrate  information  from  a variety  of 
sources.  One  example  of  die  flexibility  of  the  maps  database  is 
illustrated  by  onr  experiences  with  the  Defense  Mapping  Agency's 
(DMA)  Digital  I .andmass  Simulation  System  (DI.MS)29. 

di  ms  is  composed  of  a digiial  feature  analysis  database  (DEAD) 
which  describes  man-made  cultural  features  and  a digital  terrain 
elevation  database  (DIED)  which  is  organized  ns  a raster  elevation  grid. 
The  specified  resolution  of  the  DI-'AD  data  is  comparable  to  map  scales 
of  1:250.000  to  1:100.000.  Tile  specified  resolution  of  DIED  data  is 
within  a meter  vertical  resolution  over  a 1002  meter  (3  arc  see)  grid. 

6.4.1 . dfad:  Digital  Feature  Analysis  Database 

In  order  to  integrate  the  DI'AD  database  into  MAPS,  we  reorganized 
die  internal  DI  AD  data  structures  to  allow  for  random  access  using  a 
feature  header  list.  We  converted  die  representation  of  geodetic 
wufvtevMet,  fovm  av.  tofostk  format  that  «xx.  reLafore  to  m.  mtotwM  tat 
coordinate,  to  an  absolute  coordinate  system.  Our  dead  database 
ctv,  "rs  a two  degree  square  area,  from  latitude  N 38°  to  N 40°  and 
loiigiludt  VV  70g  ui  Vt  7$.  !i  la  Lui  jpvtSetT  of  64  roup  s’ticcis  , udl 
containing  a 1 5*x  1 5'  map  area.  We  assigned  unique  feature  identifiers 
(names)  to  map  features  because  feature  numbers  were  not  unique 
across  map  sheets.  There  arc  no  feature  names  or  semantics  associated 
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with  im:ad  entries  primarily  because  the  database  was  not  iiitcndcii  to 
be  used  as  a general  purpose  geographic  information  sy stem.  The 
feature  header  mechanism  allows  us  to  perform  random  access  to 
features  in  a map  sheet.  We  can  also  search  using  feature  attributes  such 
as  feature  analysis  code,  feature  type,  surface  material  code,  and  feature 
id  code.  This  type  of  reorganization  is  necessary  to  support  tin 
interactive  query-based  interface  for  human  and  application  programs. 

Figure  2 shows  a plot  of  polygon  features  in  the  area  corresponding 
to  our  entire  Washington  D.C.  database.  Figure  3 is  a detailed  portion 
of  the  DIM)  database  centered  on  Foggy  liottom.  For  comparison, 
l ignre  4 is  the  corresponding  area  from  the  «)\<  i i'i.mai1  database 
plotted  on  the  same  scale. 

Some  of  the  in  At)  database  entries  arc  easily  recognizable  as  natural 
or  man-made  features,  although  as  discussed,  this  information  is  not  in 
the  original  database  itself.  Figure  5 is  the  description  for  the  Tidal 
Basin,  Figure  ft  is  the  Kochambeau  Bridge.  Figure  7 is  a description  for 
a large  irregular  area  in  central  Washington  D.C.  that  contains  the 
major  government  office  holdings.  The  feature  name  assigned  by  maps 
is  tbc  first  entry  in  each  of  die  Figures. 

feature  ' <12 0 f 4 7 1 a909 ' 
feature  header:  4/1  (seek : 72416) 
feature  analysts  code:  1 002 
feature  type:  areal  feature 
surface  material  code:  (6)  water 
feature  id  code:  (909)  not  assigned 
subcategory:  fresh  water  (shallow) 
average  height  (meters):  0 

aerial  feature:  471  polygon  with  7G  vertices 
tree  cover:  0 roof  cover:  0 density:  0 
min  point  (south  west)  5298,7979 
max  point  (north  east)  5507,0385 


Figure  5:  DIM):  Description  for  Tidal  Basin 


feature  ’ d25f 4741  250 * 

feature  header:  474  (seek: 73132) 

feature  analysis  code:  1085 

feature  type:  linear  feature 

surface  material  code:  (3)  stone  / brick 

feature  id  code:  (250)  not  assigned 

subcategory ; not  assigned  (general) 

average  height  (meters):  2 

linear  feature:  474  line  with  3 vertices 

w+Hf  ■ 3 t 

first  point:  5024,8064 

last  point:  5192,8227 


Figure  6:  wad;  Description  for  Uochainbcnu  UridRC 

feature  1 d25 f 402a010 ’ 
feature  header:  402  (seek: 63688) 
feature  analysis  code:  1010 
feature  type:  areal  feature 
surface  material  code:  (3)  stone  / brick 
feature  id  code:  (01 u)  nut  uss  iyned 
subcategory  institutional  (general) 
average  height  (meters):  28 

aerial  feature:  402  polygon  with  27  vertices 
tree  cover:  10  roof  cover:  70  density:  3 
min  point  (south  west)  5705,7971 
max  poir*  (north  east)  G260.Q799 


figure  f.  OI-,vn;  Dc.Veliptofbf^OTcT™  refit  HdWbTgS 


6.4.2.  oted:  Terrain  Elevation  Database 

Tbc  organization  of  the  digital  terrain  database  is  more 
straightforward.  The  mi  l)  database  covers  tbc  same  geodetic  area  as 
our  DIM)  data.  It  is  organized  into  (>4  raster  images  using  the  same 
image  format  as  our  digital  aerial  imagery.  Each  image  containing  a 15' 
x 15'  array  of  terrain  samples,  where  each  "pixel"  is  a discrete  elevation 
point.  The  terrain  package.  Fi.PVAltOb'.  provides  a transparent 
interface  to  the  mi  l)  database.  Users  can  retrieve  elevation 
information  based  on  rectangular  geodetic  area,  closest  sample  point  to 
a geodetic  point,  or  by  weighted  interpolation.  t-t.l-VA  f ION  uses  the 
CMU  image  package  to  efficiently  buffer  blocks  of  contiguous  terrain 
data. 

6.5.  Conceptual  Map  Database 

Tbc  map  database  component  of  MAI’S,  CONCtTTMAC,  lias  been 
described  In  McKcovvir5.  We  will  give  a brief  overview  cf  the 
organization  and  concentrate  on  our  new  work  in  hierarchical 
organi/aiion  and  feature  semantics. 

6.5.1.  Concept  Schema 

The  basic  entity  in  the  conctpimau  database  is  the  concept  schema. 
The  schema  is  given  a unique  ID  by  the  database,  and  the  user  specifics 
a 'symbolic'  print  name  for  the  concept,  l-lacli  concept  may  have  one  or 
more  role  schema  associated  with  it.  Bole  schema  specify  one  ur  more 
database  views  of  the  same  geographic  concept.  For  example, 
’northwest  Washington'  can  be  viewed  as  a residential  area  as  well  as 
political  entity.  Another  aspect  is  die  ability  to  associate  the  same  name 

to  two  different  but  related  spatial  objects.  Consider  tbc  'kennedy 
center'  as  a building  ami  as  the  spatial  area  (ie.  lawn,  parking  area,  etc.) 
viiii’iiiii.iMiig  the  building,  'live  principle  nde  of  n concept  sdwms 
indicates  a preferred  or  default  view.  The  rONdl'IMAI'  database  is 
composed  of  lists  of  concept  schema. 

6.5.2.  Role  Schema 

The  role  schema  is  a further  specification  of  the  attributes  of  the  map 
feature.  It  contains  the  role  mnnr  attribute  (building,  bridge, 
commercial  area,  etc.),  a snbrale  inline  attribute  (house,  museum, 
dormitory,  etc.),  a rule  class  aunWlc  He..  tiuttaings  may  ac  gowimnna. 
residential,  commercial,  etc.),  a role  type  attribute  (ic.  physical, 
conceptual  ur  aggregate),  and  a rale  ilcrivaliun  attribute  (ie.  derivation 
method). 

Tbc  role  name,  subrule,  and  rule  class  attributes  categorize  the  map 
feature  according  to  its  function.  For  example:  this  feature  is  a 

building,  used  as  an  office  building,  used  fur  government  purposes. 
The  role  type  attribute  describes  whether  the  map  feature  is  physically 
realized  in  the  scene,  or  if  it  is  a conceptual  feature  such  as  a 
neighborhood,  political,  or  geographic  boundary.  The  role  type 
attribute  also  provides  a mechanism  to  define  the  role  schema  ns  a 
collection  of  physical  or  conceptual  map  features.  For  example,  the 
concept  schema  in  MAI'S  for  'district  of  Columbia-  has  a role  type 
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aggregraie-concepiual.  with  aggregate  roles,  ‘northwest  Washington', 
northeast  Washington',  southwest  Washington-,  and  'southeast 
Washington.  I his  mechanism  allows  the  user  to  explicitly  represent 
concepts  that  tire  strictly  composed  of  other  role  schema.  The  role 
derivation  attribute  describes  the  method  bv  which  the  role  and  its 
associated  geodetic  position  description  were  added  to  the 
roviTTMAl'  database. 

bach  role  schema  contains  a 3hi I r identifier  that  is  used  to  access  a set 
of  (ONCITIMap  database  files  which  contain  geodetic  information 
about  the  map  feature.  I hese  identifiers  can  be  shared  when  multiple 
roles  have  the  same  geodetic  description,  as  in  die  previous  example  uf 
northwest  Washington'  viewed  as  both  a residential  and  political  area. 

I he  C'ONt'LTl'MAP  313  description  allows  for  point,  line,  and  polygon 
features  as  primitives,  and  permits  the  aggregation  of  primitives  into 
more  complex  topologies,  such  as  regions  with  holes,  discommons  lines, 
and  point  lists.  Associated  with  each  feature  that  was  acquired  from  a 
image  in  tile  database  is  the  generic  name  uf  the  image.  If  the 
correspondence  of  the  generic  image  changes  due  to  the  addition  of 
more  ground  control  points,  or  better  a camera  model,  the  position  of 
the  ground  feature  can  he  automatically  recalculated. 

The  following  is  the  set  of  files  associated  with  each  linn. 

• fill]  3H  geodetic  location 

• a set  of  <latiludc/!nngitudc/clcvation>  triples  which 
deline  the  geodetic  position  of  the  role 

• [DU  | 31)  feature  shape  description 

• metric  values  for  lenglit,  width,  area,  compactness, 
centroid,  fom  icr  shape  approximation  etc. 

• (PC')  feature  image  coverage 

• a list  of  generic  images  which  contain  this  feature 

- image  nibr  and  feature  coordinates  for  each  image 

• [tutor]  feature  property  list 

• I'st  of  properties  of  the  map  feature 

- some  eenenl  (TtfUTtifl.  Uidr  **  'wjr'  Vffru'il)'  ' IH 
display  type' 

feautte  type  specific  properties  such  ,o ' uuiitbik  wf  Amiri  , 
'basement',  'height',  and  'roof  type'  for  buildings 

6.5.3.  Database  Query 

conut l map  supports  four  methods  of  database  query.  The 
methods  arc  signal  act's  symbolic  access,  template  matching  and 
geometric  access.  The  following  table  gives  a brief  description  of  each 
query  method. 

• signal  access 

Given  a geodetic  specification  (point,  line,  area)**’**, 
pvi  mi  m flit  UiHowiilg  optimums. 

• display  all  imagery  at  which  contains  point,  line  or  area. 

- retrieve  all  map  features  within  geodetic  specification 

- retrieve  terrain  elevation 


o symbolic  access 

Gi ccn  a symbolic  name,  such  as  'treasury  building'  perform 
tile  following  operations: 

- convert  inline  into  geodetic  specification  to  perform  signal 
access  operations  listed  above 

retrieve  database  description,  facts  and  properties  ol  die 
map  feature 

- retriev  e imagery  based  on  symbolic  (generic)  name 

• template  matching 

Coven  a partial  specification  of  symbolic  attributes  perform 
the  following  operations: 

- find  all  map  features  which  satisfy  the  specification 
template  and  return  their  symbolic  name 

- find  all  images  and  return  symbolic  (generic)  name 

• geometric  access 

Given  a geomctic  operation  such  as  'contains'  and  a 
geodetic  specification  perform  the  following  operations: 

1 1 nt*  all  map  features  which  satisfy  the  operation 
performed  over  the  geodetic  specification  and  return  llicir 
symbolic  name.  - find  all  image  features  and  return 
symbolic  name 

'Ihesc  primitive  access  functions  can  be  combined25  to  answer 
queries  such  as:  display  images  uf  Foggy  Bottom  before  1977',  'what  is 
the  closest  commercial  building  to  ibis  geographic  point',  and  'bow 
m;m)  nidges  cross  between  Virginia  and  the  District  of  Columbia'. 
Figure  8 is  a simple  schematic  giving  the  processes  by  which  maps 
prov  ides  signal  and  symbolic  access  into  the  concitimap  database  and 
display  of  the  query  result. 

6.5.4.  Spatial  Computation 

CONCTP'I  map  computes  geometric  properties  based  on  the  geodetic 
descriptions  associated  with  each  role  schema  in  the  database.  A static 
description  of  all  spatial  relationships  between  map  features  for 
contains,  subsumed  bv.  intersection,  ad'siccno . closest  | oint, 
partitioned  by  is  maintained  in  the  database. 

• ’contains' 

- an  unordcrcd  list  of  features  which  the  map  feature 
contains 

• 'subsumed  by 

- an  unordcrcd  list  of  features  which  contain  die  map 
feature 


Uiis  specification  may  be  in  geodclic  coordinates  or  require  imagc-lo-map 
correspondence 
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• 'intersection' 

• .in  unordcrcd  list  of  features  which  intersect  the  map 
feature 

• ‘closest  point* 

• single  feature  which  ts  closest  to  the  map  feature 

• ‘adjacent)* 

- .m  unordered  list  of  features  that  are  within  a specific 
distance  ofihc  map  feature 

• ‘partitioned  b)‘ 

• the  locus  of  points  where  two  areal  features  share  a 
common  btnindary. 

If  one  or  more  of  the  map  features  in  a spatial  compulation  is  a result  ol 
a dvnaniK  querv  (and  therefore  not  in  the  static  database),  these 
relationships  are  computed  as  needed.  N simple  memo  function  is 
implemented  to  avoid  recompilation  of  dvnamic  properties.  Ihe  use 
ol  the  static  description  can  also  he  turned  off  to  evaluate  hierarchical 
search  as  described  in  the  following  section. 


Hie  com i him  \l*  database  stores  both  factual  and  exact  inforinatiun 
describing  the  spatial  relationship.  For  example,  if  two  features 
intersect,  the  list  of  geodetic  intersection  points  is  stored,  as  well  as  the 
fact  that  they  intersect  at  least  once.  I his  is  necessary  for  query  which 
require  the  display  of  imagery  containing  a geometric  fact,  and  may 
pissihly  be  useful  for  describing  the  semantics  of  the  intersection.  In 
the  following  section  we  will  discuss  the  use  of  a hierarchical 
organi/alion  based  on  the  'contains'  relation  primitive,  and  show  how  it 
can  be  used  to  structure  the  spatial  database. 

6.6.  Fherarchical  Organization 

In  this  section  we  discuss  ihe  use  of  hierarchical  organi/alion  of 
spatial  data  in  the  M M’s  system.  Ihe  < OM  I CIM \f  database  is  used  to 
build  a hirmrcln  /nr  data  stnicturc  which  represents  the  whole-pan 
relationships  and  spatial  containment  of  map  feature  descriptions.  Ibis 
tree  is  used  in  improve  the  speed  of  spatial  compulations  by 
constraining  search  to  a pinion  of  the  database.  In  the  following 
sections  we  briefly  discuss  why  we  believe  tins  is  a govid  alternative  to 
regular  spatial  decompositions  such  as  quadtree1'  l6,  or  k-d  tree1 
usually  propised  for  MPD  model  databases. 
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6.6.1.  Regular  Decomposition 

Regular  docompusiliom  such  .is  ilu-  quadtree  organizations  Jo  not 
explicitly  exploit  the  inherent  smiclure  in  spatial  organizations. 
I’ractical  implementations  of  these  organizations  often  use  image-based 
(integer)  coordinate  systems  and  therefore  have  a hounded  position 
resolution.  In  general  cartographic  systems  it  is  important  to  he  able  to 
represent  and  manipulate  map  feature  dew  options  at  radically  diflcrent 
resolutions  using  a real  'allied  coordinate  system,  l or  example, 
consider  a dynamic  query  that  results  in  the  creation  of  a xery  small 
polygonal  area.  When  computing  containment  or  intersection  against  a 
Static  map  database  with  features  represented  as  a quadtrees,  the 
quadtrees  for  the  static  map  feature  must  be  generated  to  a much  finer 
lexcl  ofdct.nl  in  order  to  compare  the  two  data  structures.  Recent  work 
is  beginning  to  represent  quadtrees  on  real  'allied  coordinate  systems  . 
but  little  is  known  of  its  practical  implementation,  complexity,  and 
storage  efficiency.  K-d  trees  show  storage  efficiency  improvements 
over  quadtrees'7,  since  they  allow  for  a more  flexible  decomposition 
tailored  to  spatial  feature  density.  However,  they  have  tire  same 
fundarnent.il  limitations  when  used  to  represent  map  features  m a real 
valued  coordinate  swem. 

In  MAI’S  we  perform  geometric  computations  on  the  feature  data  in 
the  geodetic  coordinate  system  using  point,  line,  and  polygon  as  map 
primitives.  We  constrain  search  by  using  a hierarchical  representation 
computed  directly  from  the  underlying  map  dan.  Ilicsc  spatial 
contrumts  can  he  viewed  as  natural,  that  is.  intrinsic  to  the  data,  and 
may  have  some  analogy  to  how  humans  organize  a "map  in  the  head" 
to  avoid  search,  l or  example,  when  a tourist  who  is  looking  for  the 
Watergate  I Intel  is  told  that  the  building  is  in  Northwest  Washington, 
she  w ill  not  spend  much  time  looking  at  a map  of  \ irginia.  I depending 
on  her  familiarity  with  the  area,  she  may  avoid  looking  at  much  of  the 

map  outside  of  the  Northwest  I hsirict. As  we  begin  to  represent 

large  numbers  of  map  features  with  more  complex  interrelationships, 
we  believe  that  the  use  of  natural  hierarchies  in  urban  areas,  such  as 
political  boundaries,  neighborhoods,  commercial  and  industrial  areas, 
serve  to  constrain  search.  I hey  may  also  allow  us  to  build  systems  that 
organize  data  using  spatial  relationships  that  arc  close  to  human  spaual 

models. 

6.6.2.  Hierarchical  Decomposition 

lire  hierarchical  containment  tree  is  a tree  structure  where  nodes 
represent  map  features.  I .ich  node  has  as  its  descendants  those  features 
dial  it  completely  contains  in  clcniiuui)  space,  lire 

hierarchical  tree  is  initially  generated  by  obtaining  an  unordcred  list  of 
features  (containment  list)  for  each  map  database  feature.  Starting  w ith 
a designated  root  node  ('greater  Washington  d.c.  ) which  contains  all 
features  in  lire  database,  descendant  nodes  are  recursively  removed 


******  If  eic  tv  told  lhai  Ihc  Watergate  tv  also  near  die  Potomac  nver.  that  should 
further  constrain  her  scaich.  bill  lhal  is  another  story 


from  the  parent  node  list  if  they  are  already  contained  in  another 
descendant  node.  The  result  is  that  the  parent  node  is  left  with  a list  of 
descendant  features  that  are  not  contained  by  any  other  node.  Ihcsc 
descendant  nodes  form  die  next  level  of  an  N-.iry  tree  ordered  by  the 
contains  relationship.  This  procedure  is  performed  recursively  for 
every  map  feature.  Terminal  nodes  arc  point  and  line  features,  or  areal 
features  that  contain  no  other  map  feature.  We  will  discuss  the  point 
containment  and  closest  point  compulation  using  the  hierarchy  tree  in 
the  follow  ing  section. 

figure  9 shows  a small  section  of  the  hierarchical  containment  tree, 
lire  use  of  conceptual  features  - features  with  no  physical  realization  in 
the  world  but  represent  well  understood  spatial  areas-  can  be  used  to 
partition  the  database.  In  this  case  the  map  feature  foggy  bottom' 


4!  entries  lor  ‘contains'  for  ‘northwest  Washington’ 
entry  O'  ‘ncnillan  reservoir  (role:  0)’ 
entry  1 kennedy  center  (role:  0)' 

entry  2 ’ellipse  (role  0)’ 

entry  3 executive  office  building  (role  0)‘ 

entry  4 white  house  (role-  0)‘ 

entry  b:  'treasury  building  (role:  0)' 

entry  f.  ‘department  of  commerce  (role  0)' 

entry  / museum  of  Instory  and  technology  (role:  0)' 
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'constitution  avenue  (role 
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'national  cathedral 
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entry 
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40 
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or 

\\  entries  for  ’contains'  for  'foggy  bottom' 
entry  0 'ker  edy  center  (role:  1)' 
entry  1 Washington  circle  (role  0)' 
entry  2 slate  department  (role.  0)’ 

entry  3:  amen  can  pharmaceutical  association  (role:  0) 

entry  4 national  academy  of  sciences  (role  0)' 

entry  b federal  reserve  hoard  (role  0)* 

entry  f>  national  science  foundation  (role:  0)' 

entry  7:  'civil  service  commission  (role  0)* 

entry  8 ‘c  street  (role  0)' 

entry  1:  '??nd  street  (role:  0)’ 

entry  10  'south  new  hampshire  avenue  (role  0)* 


Figure  9:  MAPS:  Hierarchical  Spatial  Containment 
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allows  us  lo  partition  some  of  the  buildings  and  roads  dial  are  contained 
within  'northwest  Washington',  As  more  neighborhood  areas  and  city 
districts  arc  added  to  our  database,  we  expect  to  see  improved 
performance  especially  in  areas  with  dense  feature  distributions.  This 
will  also  improve  the  richness  of  the  spatial  description  available  to  the 
user. 

6.6.3.  Hierarchical  Search 

In  this  section  we  discuss  the  use  of  our  hierarchical  organization  to 
partition  the  map  database  to  improve  performance  by  decreasing 
search  when  computing  the  spatial  relationships  of  map  features.  Ibc 
hierarchical  searching  algorithm  is  basically  an  N-ary  tree  searching 
algorithm.  Consider  a user  at  the  concii'IMAI’  image  display  who 
invokes  die  gcometic  database  to  compute  a symbolic  description  of 
what  map  feature  he  is  pointing  at.  First,  using  image-to-map 
correspondence,  die  system  calculates  the  following  map  coordinates: 

latitude  N 38  53  40  (276) 

longitude  W 77  03  53  (337) 

This  point  is  converted  into  a temporary  map  database  feature  and  is 
tested  against  the  root  node  of  the  hierarchy  ire,.  If  it  is  not  contained 
in  this  node  (not  generally  the  case),  men  toe  point  cannot  correspond 
to  a database  feature,  and  the  search  terminates.  The  user  is  informed 
that  the  point  is  outside  the  map  database.'*""”  If  the  'contains*  test 
succeeds,  it  rectirses  down  the  tree  and  performs  the  test  against  the 
siblings  of  the  node  just  tested.  The  search  allows  several  paths  to  exist 
for  any  point,  dins  more  than  one  sibling  may  contain  a path  to  the 

point.  Ibis  soil  of  aiioiiial)  oceuis  vt1u.ii  a Kaltiic  to  u-Mat  ill 

the  intersecting  region  of  two  larger  regions.  However,  if  the  feature  is 
not  contained  by  the  node,  it  is  not  contained  by  any  of  the  node's 
Ueseerrfsrttts-.  itftd  that  portion  of  the  tree  ir  ir-tt  ftttflw  setwe?texf.  Fign'ti 
10  shows  the  answer  to  our  hypothetical  query.  The  query  point  is 
contained  wi'hin  Theodore  rooscvclt  island',  and  two  search  paths  in 
the  owtainmeoi  tree  arc  given  The  same  mo+vanism  is  used  ft'1- lin'' 
and  polygon  features,  although  the  primitive  determination  of 
containment  depends  on  the  geometric  type  of  the  feature. 


6.7.  Toward  Feature  Semantics 

We  have  begun  to  investigate  the  generation  of  map  feature 
semantics  directly  from  die  hierarchical  representation  of  the  map 
feature  data.  A simple  example  is  the  semantic  description  of  a bridge: 
the  feature  names  and  map  locations  that  it  connects  as  well  as  die 
names  of  the  map  features  diat  it  crosses  over.  Figures  1 1 and  12  show 
the  result  of  applying  a procedural  description  of  the  semantics  of  a 
bridge  concept  to  calculate  the  'connects'  and  'crossover'  relationship 
using  the  map  feature  descriptions  of  'aldington  memorial  bridge'  and 
Theodore  rooscvclt  memorial  bridge'.  These  results  are  generated 
directly  using  the  M \l>s  hierarchical  organization  for  spatial  data.  We 
do  not  pose  this  as  a theory  of  map  feature  semantics,  but  envision  a set 
of  feature  specific  procedures  that  can  build  these  types  of  descriptions. 


2 entries 

f or 

’contains'  for  'querypoint  1' 

entry  0: 

'Virginia' 

entry  1: 

'greater  Washington  d.c.’ 

* 

2 entries 

for 

'contains'  for  'query point  1' 

entry  0: 

'arlington  memorial  bridge 

entry  1: 

'greater  Washington  d.c.' 

entries  for 

contains  or  querypoinl  2 

entry  0: 

'mall  area' 

entry  1 : 

'southwest.  Washington' 

entry  2: 

'district  of  Columbia’ 

entry  3; 

'greater  Washington  d.c.’ 

2 entries  for 

'contains'  for  'querypoint  2' 

entry  0: 

'arlington  memorial  bridge 

entry  1: 

'greater  Washington  d.c.’ 

5 entr 

ies  for 

'intersection'  for  'crossover 

entry 

0: 

'Virginia* 

entry 

1: 

'district  of  Columbia' 

entry 

2: 

'southwest  Washington' 

#*lFf 

Jii 

’Ml'  r.  * 

entry 

4; 

'potomac  river  {Role:  0)' 

2 entries  fur  'connects'  for  'Arlington  memorial  bridge 
entry  O'  ' v i r ( in i a ' 

entry  1:  'mall  area' 

1 entries  for  'crossover'  for  'arlington  memorial  bridge 
entry  0:  'potomne  river’ 


This  node  belongs  in  the  following  place(s): 

3 entries  for  'contains'  for  ‘theodore  roosevelt  island' 
entry  0:  'northwest  Washington’ 

entry  1:  'district  of  ‘Columbia' 

entry  2:  'greater  Washington  d.c.’ 

AND  •• 

2 entries  for  'contains'  for  'theodore  roosevelt  island' 

entry  0:  'potomac  river’ 

entry  1:  ’greater  Washington  d.c.' 


Figure  10:  MAI'S:  Containment  T ree  linlry  fur 
T heodore  Roosevelt  Island 


flfis  can  actually  octui  i-  c.ira  ■X’CTTnr.  ft/tmTcrtVS 

trough  i he  terminal.  Therefore  the  dalabase  has  some  crude  idea  of  its  cxlcnl  of  map 
knowledge 


Figure  II:  maps:  Semantic  Computation  from  Spatial  Data 
Arlington  Memorial  Bridge 

T he  procedure  for  bridge  semantics  is  as  follows:  A bridge  can  be 
represented  in  die  CONCTT’TMAI'  database  ns  an  polygonal  area,  a list  of 
linear  segments,  or  as  a geodetic  point.  The  polygonal  area  arises  when 
die  bridge  deck  is  represented,  die  list  of  linear  segments  approximates 
die  center  line  of  the  bridge,  and  the  point  feature  generally  represents 
that  die  bridge  is  a landmark  feature.  No  semantics  are  computed  in 
the  latter  case.  If  the  bridge  is  represented  as  a line,  the  end  points  are 
selected,  otherwise  the  endpoints  of  the  major  axis  of  the  bounding 
ellipse  arc  -ctrievcd  from  the  feature  ()31T  file.  At  some  level  of 
description,  diese  endpoints  define  the  ‘connects*  relationship,  but  this 
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2 entries  for  'contnins*  for  ’i|iierypoint  1* 

entry  0:  ’Virginia’ 

entry  is  ‘greater  Washington  d.c.* 

.............  AND  MMM.MMM 

2 entries  for  ’contains'  for  'querypoint  1’ 

entry  0:  ’theodore  roosevelt  memorial  bridge' 

entry  1:  ’greater  Washington  d.c.’ 


3 entries  for  ’contains*  for  ’nuerypoint  2* 
entry  0:  ’northwest  Washington’ 

entry  1:  'district  of  Columbia* 

entry  2:  greater  Washington  d.c.’ 

A N 0 •••••*••*•••• 

2 entries  for  ’contains*  for  ‘querypoint  2* 

entry  0:  'theodore  roosevelt  memorial  bridge' 

entry  1:  'greater  Washington  d.c.’ 


£>  entries  f nr 

'intersection*  for  'crossover  list* 

entry  0: 

’Virginia’ 

entry  1; 

’district  of  Columbia* 

entry  2 : 

'northwest  Washington’ 

entry  3: 

’theodore  roosevelt  island’ 

entry  4: 

'potomac  river' 

2 entries  for  ’connects'  for  'theodore  roosevelt  memoria 

entry  0:  'Virginia* 

entry  1;  'northwest  Washington* 

2 entries  for  'crossover'  for  ’theodore  roosevelt  memori 
entry  0:  'theodore  roosevelt  island* 

entry  1:  'potomac  river* 


I'ignrc  12:  maps:  Semantic  Computation  from  Spatial  Data 
Theodore  Koose\clt  Memorial  Bridge 
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is  not  useful  if  we  arc  envisioning  generation  of  a reasonably  complex 
symbolic  representation. 

Tlie  'contains'  relationship  is  applied  to  each  endpoint  using  the 
hierarchical  tree  to  order  die  search.  As  before,  this  search  returns  a list 
of  features  ordered  by  spatial  eoM.ihvmonl,  and  there  may  be  several 
independent  containment  paths.  Redundant  paths  arc  eliminated  by 
examining  whether  die  bridge  is  in  die  containment  path.  The  first 
entry  (It)  in  each  of  the  remaining  paths  is  one  of  the  areas  connected  by 
the  laidge.  Using  the  'contains'  relationship,  the  other  entries  in  the 
path  are  also  valid  connecting  areas. 


7.  Synthesis  Tasks 

In  this  section  we  will  discuss  three  applications  of  the  MM'S  database 
to  cartographic  and  image  interpretation  tasks.  These  tasks  arc  3D  scene 
generation  of  views  of  Washington  I).  C„  die  use  of  the  map  database 
to  guide  image  segmentation,  and  some  preliminary  results  on  a rule- 
based  sjmciii  for  airpuii  scene  interpretation.  Ivacfi  task  requires  the 
capabilities  of  various  aspects  of  the  IMIi  model  as  implemented  in  the 
MAI'S  system.  These  applications  pull  together  external  and  image/map 
latabascs.  and  are  only  possible  using  an  integrated  system  th  it  relates 
imagery,  terrain,  and  map  data  dirough  a unified  cartographic 
representation. 


To  compute  die  'crossover'  relationship,  the  'intersection' 
relationship  is  computed  for  the  bridge  using  die  complete  list  of  line 
segments  or  the  polygonal  description.  A list  of  all  the  features  that  die 
bridge  intersects  is  assembled.  1-ntries  in  the  intersection  list  are 
removed  if  tiicy  arc  alio  prcscni  in  either  oi  me  eonfieeiS  lisi.v  Tne 
assumption  is  that  those  features  that  didn't  contain  a bridge  endpoint, 
hut  in lorsccted  with  the  bridge  description,  arc  those  features  that  the 

hudi,c  dosses  ovo.  II  dieie  11  sufficiently  detailed  cievudon  data  fui 
man-made  features  it  should  be  possible  to  compute  semantics  for 
p 1-  ■ 'r  iw-tr'  p.vtv«  iiislw'  fiy  eliiinihiT-a-g  ■he  h*  -it  st 

die  actual 


7.0.1 . WASH3D:  3D  Scene  Generation 
Hie  first  application  of  the  MAI'S  database  is  in  the  area  of  315 
computer  graphics  for  scene  simulation  and  database  validation. 
Computer  graphics  play  an  important  role  in  the  areas  of  image 
piwicssmg.  ptioui  micijucuiiimi,  unu  cauogi aptly.  In  eu. tugrapiiy 
various  phases  of  the  map  generation  process  use  graphics  techniques 
or  source  matciial  analysis,  transcription  and  update,  and  some  aspects 
of  nmp  myuut  mid  pioJutAiou.  lluwcvtf.  Ilumy  iilajol  amps  in  Hie 
generation  of  a cartographic  product  remain  largely  manual.  One 
'm’purriw  fcf  ishiiJh  hmlev/nifte  tools  edit  is  the  integral..  ■.  .1 

is  often 


Figure  14:  WASII3D:  Vertical  View  85°  Northwest  Washington 
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• igurc  15:  WASII3D:  Northwest  Washington  Front  Above  National  Airport 


usetl  to  verily  the  geodetic  accuracy  of  natural  and  man-made  features 
in  tile  digital  database  prior  to  actual  map  layout  .md  production. 
Another  application  is  sensor  simulation’0  1I,  Radar,  v isual,  and  multi- 
sensor  scenes  arc  digitally  generated  to  verify  the  quality  of  digital 
curtoic  and  lenain  databases  oi  to  dcicoimic  die  quality  of  tile  sensoi 
model.  Improvements  to  the  level  of  detail  contained  in  the  underlying 
database  can  be  subjectively  measured  in  terms  of  the  quality  of  the 
getter, tied  scene. 

tt'  i 'iih"1’*  rwMHwniwhggM*-.  flt-iystertMto*  n*.r  the  M.M 
to  integrate  a digital  terrain  database,  a cultural  feature  database,  and 
the  cost ’t  ft \t.M'  database  to  allow  a user  to  generate  cartographically 
accurate  rues  for  human  visual  analysis,  wash  in  uses  the  coarse 
resolution  DIMS  database  described  in  Section  6.4  to  generate  a 
h tselmc  thematic  map.  The  thematic  map  is  a 21)  image  which  is 
produced  by  scan  conversion  of  the  ni  MS  digital  feature  analysis 
lUMhaa*  <f»  .puf?#uu  J*U  UlJi*  We  j*>iijH  s Wtltw  un  rsgiw 


water  (blue),  residential  (yellow),  and  high-density  urban  (brown), 
in  ms  terrain  elevation  dam  (mm)  is  interpolated  to  determine  ground 
elevations  at  each  point  in  the  21)  image.  Since  the  resolution  of  die 
in  ad  data  is  coarse,  comparable  to  map  scales  of  1,250,000  to 
1: 1IRJ.UUU.  vve  use  die  i u i i i.vmr  dat.itrase  to  pfov.de  high  TCsakmor. 
31)  feature  descriptions  of  buildings,  roads,  bridges,  residential  and 
commercial  areas.  The  cosci  I'lMAI’  database  is  derived  from  imagery 
with  resolutions  between  1:12000  and  1:36000.  and  die  addition  of 
these  features  effectively  intensifies  die  perceived  level  of  detail  in  die 

siimu.uvu  scene,  e en  a >wLi.a i a iv  k-,...c  tn'i  bv  lit 


II  {JlUwiJl  ('lie  br,.:c  rr„  , 


I . tikes'1  describes  the  utility  of  selective  database  intensification  for 
tailoring  standard  database  products  to  custom  applications  and  for 
tiinc-ci ideal  applications  which  cannot  be  handled  by  normtil 
pioduciion  schedules,  Figure  13  ihows  the  imcractlve  pweCTS  by  which 
users  can  specify  an  area  of  interest  for  31)  scene  generation.  Figures 

Ttti  -l!  tlhrW'  hwi  Mu'  Wwdblig I 11  alerl  lu'immied 

hv  WASII3D. 
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7.0.2,  machineseg:  Map-Guided  Machine  Segmentation 
The  second  application  of  the  maps  datahase  is  in  the  area  of  map- 
guided  machine  sOgiiieiuauoii.  fj.ssis  may  ''pL,.d\  .*  I II..  I’  fc*.Vuu-  if  on  I 
the  tOsiTI’IMAl*  database  or  interactively  generate  a feature 
vliM-tlpHtth  a sing  *c  r program.  fcr  ww-hm*  w*  map  dws*»e 

feature.  MACIllM  Sl  fi  uses  an  existing  image  coverage  |I  C|  tile  (sec 
SojtiuH  h.S  ?)  that  ipixUkt  hi  wliirh  imfes  Hie  feature  is  found,  and 
the  feature  location  in  the  image,  l-'or  interactive  specification,  an  |l  C) 
file  is  created  dynamically  by  image-lo-map  correspondence  using  the 


l-or  each  image,  a high  resolution  window  containing  die  database 
ic.nme  n , Al I .Il  ls d and  displayed.  We  eepand  die  sl/.e  of  die  linage 
window  to  contain  an  area  of  uncertainty  around  the  feature  location. 
'Ifw  c vpiiii.tiun  l!i  aitHMiitk  ksirv!  yjrt  iW vbe  ftWVffC.  H«  -we  fmtl 

to  incorporate  correspondence  error  measures  based  on  the  quality  of 
'ha  ('inner, i model  associated  with  each  image  The  image  window  is 
smoothed,  and  a segmentation  is  performed  using  a region-growing 
technique1'1  which  combines  an  edge  sir  ngth  metric  and  region  merge 

uctepttW'r.  based  an  npeetefl  wnrllfirte  !o  CuiWflJ  teghm  gmwrng, 


l-'igurc  16  shows  die  segmentation  of  several  low-elevation  buildings 
along  the  perimeter  of  die  Washington  Kllipsc.  The  uppermost 
building  is  added  to  the  CO\rmMAP  dalnhasc  in  the  standard  manner 
described  in  Section  6.5.  The  user  specifies  the  image,  iXJ8tii7.  to 
perform  the  ’^mentation  and  the  mauumsi  o system  .luuim.iiically 
displays  a reduced  resolution  window  of  the  image  Ulc3SHI7),  and  a 
high  resolution  window  ( ellipse  urea)  containing  the  database  area. 
MAf  lltM  si  c creates  a copy  of  the  high  resolution  window  as  a work 
area  I \n  ,ni,lc | for  fix  idwjHr  PitceeMfuJ  nomine#  "0  JfU-'4§l'  «iniHirfiihc 

operation  is  followed  by  the  generation  of  seed  regions  using  a 
Conservative  similarity  measure  in  insure  that  potentially  matchablr 
regions  arc  not  prematurely  merged.  The  initial  seed  regions  arc 
overlaid  on  the  image  using  graphics  overlays.  Any  seed  regions  that 
satisfy  the  shape  criteria  for  die  databac.  feature  arc  extracted  and 
marked.  In  this  example,  the  database  feature  itself  was  marked  in  die 
initial  seed  region  matching.  As  regions  arc  merged  based  on  weak 
edg^  i ouridtineS  atiJ  Tiigh  .spectral  eotiipiiuTiJTiiy,  die  iCSuTdog  icgion  is 
evaluated  with  respect  to  a list  of  shape  and  spectral  criteria.  If  the 
region  satisfies  the  criteria,  it  is  marked,  and  further  merging  is  allowed 
only  if  the  proposed  merge  improves  the  overall  region  score.  Criteria 
include  fractional  fill.  area,  linearity,  perimeter,  compactness,  and 
spectral  measures. 

The  final  results  arc  shown  in  the  second  window  labeled  set  aside. 
five  buildings  similar  to  the  map  database  feature  were  correctly 
identified  while  one  building  was  omitted.  Six  segments  were 

incorrectly  identified.  I lad  we  made  use  of  spectral  information  in  this 
particular  segmentation-  that  the  building  roofs  were  bright 
features-  we  probably  could  have  excluded  5 of  the  6 errors.  However, 
we  are  more  concerned  with  using  weak  knowledge,  and  one  cannot 
expect  better  performance  without  more  sophisticated  analysis. 
MACIIIM:si:g  allows  the  user  to  delete  erroneous  segments  and 
generates  map  descriptions  of  each  extracted  feature.  These 
descriptions  can  then  he  used  to  search  for  these  features  in  other 
database  imagery. 

The  significance  of  maciiistspu  is  that  it  can  search  systematically 
for  features  in  a database  of  images,  an  operation  that  is  fundame  ,ul 
for  change  detection  applications.  It  djrcclly  uses  tile  map  database 
description  as  an  evaluation  tool  for  image  segmentation  and 
interpretation.  It  also  uses  very  general  image  processing  tools  to 
perform  both  segmentation  and  evaluation  and  is  amenable  to 
supporting  other  approaches  to  image  segmentation  and  feature 
recovery.  A further  application  of  the  maciiinpspg  system  is  discussed 
in  the  following  section. 

7.0.3.  sr  am:  Rule-based  System  for  Airport  Interpretation 

The  third  application  of  the  MAI’S  system  is  in  the  investigation  of 
rule-based  systems  for  the  control  of  image  processing  and 
fflicn-Tsteh  wtSi  respect  t>,  s vrarfd  ti  trdr1. 


In  photo-interpretation,  knowledge  can  range  from  stereotypical 
information  about  man-made  and  natural  features  found  in  various 
situations  (airports,  manufacturing,  industrial  installations,  power  plnnts 
etc.)  to  particular  instantiations  of  these  situations  in  frequently 
tViuiiitoied  sites.  It  is  crucial  Tin  photo-interpretation  applications  that 
the  metrics  used  be  defined  in  a cartographic  coordinate  system,  such  as 
<huiiutlc/hmgiiude/cle\aiiun>.  rather  than  an  image-based  coordinate 
system.  Descriptions  such  as  "the  runway  has  area  12000  pixels"  or 
-m’  !#.*«t  'H'»  and  J4j>  p4«k"  are-  wietes  fr 
(perhaps)  the  analysis  of  one  image.  It  is  the  ease,  however,  dial  to 
operationalize  metric  knowledge  one  must  relate  the  world  model  to  the 
image  under  analysis.  This  should  lie  done  through  imngc-to-map 
correspondence  using  camera  models  which  is  the  method  used  in  our 
system. 

We  have  begun  to  build  spam15  to  test  our  ideas  in  die  use  of  the 
cu/TVb.iiiiTHiri  ii  .i  nuip  uiiiiib.se,  task  aulTe pel ijcik  low-level  image 
processing  tools,  and  a rule-based  system. 

spam  uses  the  maps  database  to  store  facts  about  man-made  or 
natural  feature  existence  and  location,  and  to  perforin  geometric 
computation  in  map  space  rather  than  image  space.  Differences  in  scale, 
orientation,  and  viewpoint  can  be  handled  in  a consistent  manner  using 
a simple  camera  model.  The  maps  database  facility  also  maintains  a 
partial  model  of  interpretation,  separate  from,  but  in  die  same 
representation  a's.  the  map  feature  database. 

The  image  processing  component  is  based  on  die  maciiintsi-G 

ptMnT MV,  uOSeTibcJ  tVi  the  pluvious  section.  It  p'crfoiTns  low  level  and 
intermediate  level  feature  extraction.  Processing  primitives  arc  based 
on  linear  feature  extraction  and  region  extraction  using  edge-based  and 
region-growing  techniques.  It  identifies  islands  of  interest  and  extends 
those  islands  constrained  by  die  geometric  model  provided  by  maps 
and  model-based  goals  established  by  die  rule-based  component. 

I lie  rule-based  component  prov  ides  the  image  processing  system 
with  the  best  next  task  based  on  die  strcngdi/proniisc  of  expectations 
and  with  constraints  from  die  image/map  database  system.  It  also 
guides  the  scene  interpretation  by  generating  successively  inure  specific 
expectations  based  on  image  processing  results. 

We  arc  in  tile  preliminary  stages  of  development  for  die  spam  system 
and  have  begun  to  build  a detailed  map  model  of  National  Airport, 
l-'igurc  17  gives  an  example  of  die  ability  of  the  maps  database  to  use 
image-to-map  correspondence  to  generate  unified  spatial  models  from 
partial  information.  The  line  drawing  labeled  J740I.IMC  contains  the 
northern  section  of  National  Airport;  M09IMG  is  a partially 
overlapping  southern  section  of  National  Airport.  Line  segments 
represent  point,  line,  and  areal  features  corresponding  to  runways, 
terminal  buildings,  access  roads,  and  hangars,  interactively  specified 
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Kigiirc  17:  spam:  Naiiomil  Airport  Spatial  Model 


unilied.img 


Unified  scene  ol  Washington  National  Airport 


using  the  CON'dil’lMAl’  rcpresenlalion.  For  those  features  that  appear 
i't  hotlt  images,  the  concept  rob  mechanism  (see  Section  6.5,2)  is  used 
to  specify  multiple  (Inlilmlc/longilittb/cbnilion}  descriptions.  A 
unified  map  description  is  created  by  matching  corresponding  line 
segments  using  the  overlapping  image  areas  (in  map  space)  to  constrain 
search,  The  result  of  unification  is  die  line  drawing  hibclcd 
AlRI’ORT.lMG. 

8.  Future  Work 

Our  future  work  will  be  directed  toward  two  research  topics.  First, 
we  have  only  begun  to  explore  die  use  of  MAI’S  as  a component  of  an 
image  interpretation  system.  We  will  continue  our  work  in  die  airport 
scene  interpretation  task,  using  the  spam  system  as  a testbed  Tor 
integration  of  a rule-based  system  with  the  maps  system.  Second,  there 
is  much  to  do  in  expanding  die  CONCFI’IMAp  database  to  include  more 
complex  3D  descriptions,  and  in  attendant  issues  or  scaling  and  si/.ing 
to  larger  databases.  Other  tasks  we  will  pursue  arc  the  evaluation  of  our 


Northern  section  ol  Washington  National  Airport 


Southern  section  ol  Washington  National  Airport 

hierarchical  spatial  representation  to  constrain  search  in  large  databases, 
general  solutions  to  complex  spatial  queries  for  situation  assessment 
applications,  and  the  application  of  spatial  knowledge  to  navigate 
through  a map  database. 

In  discussing  future  work  it  is  important  to  understand  the  strengths 
mid  limitations  of  the  current  research.  Tile  strengths  of  this  work  lie  in 
several  unique  features  of  die  maps  system.  First,  we  have  constructed 
a system  of  moderate  complexity  which  lias  signincam  capabilites  in 
each  area  of  our  Image/Map  Database  model.  The  system  integrates 
map  knowledge  from  diverse  sources  and  performs  several  tasks  diat 
require  synthesis  of  dus  knowledge.  We  have  die  ability  to  represent 
complex  map  features  in  a uniform  cartographic  coordinate  system  and 
can  compute  new  spatial  relationships  directly  from  the  map  data. 
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The  major  limitation  in  the  maps  system  is  the  current  method  for 
performing  image-to-ntap  correspondence.  from  the  standpoint 

of  the  state  of  the  art  in  photogrammetry.  we  make  simplistic 
plancnictric  itssttmptions  in  our  correspondence  algorithm,  hut  they  do 
give  reasonable  results  for  several  reasons,  Hirst.  all  of  our  photographs 
are  vertical  aerial  mapping  imagery,  and  eflorts  are  taken  to  minimize 
camera  tilt.  Second,  we  have  very  high  resolution  photographs,  each  of 
which  covers  a relatively  small  atm  ami  due  to  the  relatively  local  level 
terrain  in  Washington  I).  C„  our  polynomial  correspondence  functions 
arc  reasonably  accurate. 

The  issue  is  not  how  to  recover  camera  information  from  the 
imagery,  since  in  cartography  and  manual  photo-interpretation  the 
sensor  models  and  ephemeral  data  arc  vvcll  known  and  modeled,  but  to 
use  existing  photograminctric  tools  for  basic  data  acquisition. 
Therefore,  in  this  limitation  we  sec  an  opportunity  to  investigate  how 
MAPS  could  he  interfaced  to  a photograminctric  frontend  which  would 
Jhccitv  provide  ffr/tfff.dl^/rTgW.WvfP^1 1 Afltt  tow*  A 
model."”"'"  The  frontend  should  have  a landmark  database  and 
interactive  display  tools  to  guide  the  stereo  model  setup  in  a manner 
simiLr  v.  out  ctmem  implement  .stiwt.  Nothing  in  the  current  maps 
implementation  precludes  such  an  interface  since  we  maintain  a 30 
map  feature  representation  throughout  tire  database  using  the  USGS 
terrain  database.  The  building  of  such  tools  should  lie  the  common 
objective  both  to  cartographers  and  to  computer  scientists. 
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maps  System  Major  Components 


This  Appendix  eonlams  a list  of  Hie  major  program  modules  which 
compose  lire  maps  system. 


NAME  Size  (bytes)  COMMENTS 


Browse 

browse 

305500 

interactive  image  display  facility 

picpac 

530762 

interactive  image  processing  facility 

Corres 
cor  res 

206042 

interactive  image-map  correspondence 

checkcorres 

49523 

check  correspondence  errors 

cormain 

52093 

correspondence  algorithm 

corpa irs 

751 76 

edit  correspondence  pairs  file 

creatsdf 

50601 

create  a scene  description  Tile 

duinpeoef 

19649 

dump  a coefficients  file 

dumpeor 

23547 

dump  a correspondence  file 

dumpsdf 

25390 

dump  a scene  description  file 

hypeorpa irs 

023B0 

generate  hypothesized  landmarks 

updatesdf 

59099 

update  a scene  description  file 

Landmark 

1 nmlmark 

194953 

interactive  landmark  extraction 

creat 1 rim 

23557 

create  binary  landmark  file 

uty tod3 

5021  7 

make  a .03  file  from  an  .ety  Tile 

oty tol dm 

19948 

create  landmark  file  from  .ety  files 

1 (laser  ibe 

43695 

give  landmark  descriptions 

1 dinr  i prt 

30275 

dump  all  info  about  a landmark 

ldmtes  t 

28696 

find  landmarks  within  geodetic  area 

Seyment 
segment 
ink  idf 

170230 

hand  segmentation  program 

10537 

create  ascii  Tile  from  binary  seg  fil( 

seg rename 

39045 

edit  segmentation  region  names 

Machineseg 

nar.  h \ nosey 

290222 

machine  segmentation  program 

Conceptmap 

conceptmap 

665710 

associate  conceptual  and  map  data 

hu  i ldsegntap 

90301 

build  composite  segmentations 

coe  track 

125241 

track  points  using  map  correspondence 

congeoal 1 

21327B 

generate  geometric  database 

d 3 dump 

24629 

dump  o 03  file 

d3entcor 

93936 

create  corres  entry  from  .d3  file 

d3f  dump 

31039 

dump  a d3  f eature  file 

d3tod3f 

15826 

convert  a .03  file  to  a feature  file 

d3  to  img 

44710 

generate  binary  image  from  . d3  files 

dlmsseg 

12B324 

create  DIMS  overlay  for  geodetic  area 

dmaex tract 

31544 

extract  features  from  DLMS  . fea  files 

dunpql 

207962 

dump  a query  1 ist  file 

dumpsdf 

25390 

dump  a scene  description  file 

ecdump 

9425 

dump  the  contents  of  a coverage  file 

ecshow 

137700 

display  manager  for  coverage  files 

ecsor t 

26624 

sort  coverage  files  by  keys 

ectcseg 

18173 

create  .seg  file  from  coverage  file 

hierarchy 

486262 

huild  and  access  hierarchical  database 

h iertrack 

321869 

track  and  display  pts  using  hierarchy 

idh ier 

254739 

identify  points  using  hierarchy 

imagetoec 

34203 

associate  image  with  coverage  file 

imagetomap 

54092 

<gener  ic>«.rowxcol>  °>  <1  at/ lon/el ev> 

photo 

299710 

interactive  image  pho togrammetry 

segtod3 

57034 

convert  .seg  file  to  .d3  data  structure 

sogtoimg 

32785 

convert  .seg  regions  to  binary  image 

stereoshow 

153125 

show  stereo  image  pairs 

un i fyseg 

107603 

unify  segmentation  regions 

Wash3d 

wash3d 

764517 

3d  scene  generation  from  MAPS  database 

df  eapr t 

45013 

print  DLMS  feature  given  dims  code 

d ispfea 

137335 

display  a DLMS  map  feature  file 

dims 

53134 

create  dims  index  file 

dlmsbin 

34604 

convert  ascii  feature  files  to  binary 

dlmsf ind 

45419 

find  a DLMS  feature  based  on  attributes 

f eadumper 

45267 

dump  a 0LMS  feature  file 

Terrain 

elevation 

24097 

access  terrain  data  images 

127 


■ a&a.. 


